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A coherent state representation for the electrons of ordered antiferromagnets is used to derive 
effective Hamiltonians for the dynamics of holes in such systems. By an appropriate choice of 
these states, the constraint of forbidden double occupancy can be implemented rigorously. Using 
these coherent states, one arrives at a path integral representation of the partition function of the 
systems, from which the effective Hamiltonians can be read off. We apply this method to the t-J 
model on the square lattice and on the triangular lattice. In the former case, we reproduce the 
well-known fermion-boson Hamiltonian for a hole in a collinear antiferromagnet. We demonstrate 
that our method also works for non-collinear antiferromagnets by calculating the spectrum of a hole 
in the triangular antiferromagnet in the self-consistent Born approximation and by comparing it 
with numerically exact results. 



I. INTRODUCTION 



It is widely recognized that an accurate description of 
the dynamics of holes in antiferromagnetically (AF) or- 
dered materials is an important first step towards an un- 
derstanding of the essential physics of the cuprate suaec^ 
conductors. Consequently, numerous investigationsaO 
have dealt with the problem of a single hole in a square- 
lattice antiferromagnet. Most ofJicse studies are based 
on an effective Hamiltonian 7i n ll3llil which describes the 
hole as a spinless fermion which is coupled to bosons rep- 
resenting the collective excitations (magnons) of the anti- 
ferromagnetic background. The structure of this effective 
Hamiltonian reflects the two-sublattice structure of the 
square lattice AF. When the hole hops between nearest 
neighbour lattice sites, it necessarily disturbs the mag- 
netic order. Thus, any move of the hole requires absorp- 
tion or emission of a magnon. Derivations of Tin from the 
t-J Hamiltonian have revealed the approximate nature 
of this effective Hamiltonianlij, and to obtain the spec- 
trum of a hole from TL a further approximations are neces- 
sary. Most often, the fermion-magnon coupling is treated 
in the self-consistent Born approximation (SCBA) which, 
was already used in tie original studies by Kane et al.E3 
and Marsiglio et al.E-3. Surprisingly, the approximate 
spectra obtained in this manner agree excellently with 
the spectra of small clusters obtained by numerically ex- 
act diagonalization of the t — J Hamiltonian for these 
clusters. 

Apart from its relevance to the physics of cuprate su- 
perconductors, the dynamics of holes in ordered antifer- 
romagnets is a highly nontrivial problem of the physics 
of strongly correlated electron systems, and it is there- 
fore of fundamental theoretical interest. Here, we shall 
present a derivation of effective Hamiltonians for holes 



in AFs which is not restricted to any particular type of 
magnetic order. Our method does not provide us with a 
unique effective Hamiltonian for an arbitrary number of 
holes. Rather the cases of a single hole, two holes and in 
general an arbitrary but fixed number of holes have to be 
treated separately. While our presentation remains gen- 
eral, we shall confine ourselves to the single-hole problem 
in order to demonstrate our method's practical applica- 
bility. In particular, we shall derive and evaluate the 
effective Hamiltonian Ti a for a single hole in the trian- 
gular AF. 

II. EFFECTIVE LAGRANGIAN 

For definitcncss, we consider the t — J Hamiltonian, 
H = -tJ2 P{cL a c r ,. a + h.c.)P 

<r,r'> 

a 

+ J Sr-S r ,. (1) 

<r,r'> 

t and J arc the single particle hopping matrix clement 
and the exchange constant, respectively, the operator P 
projects onto states in which each of the N s lattice sites is 
either empty or singly occupied, and S r = \c\ a cr a ^c r a 
is the electron spin expressed in terms of creation and 
annihilation operators, cj a and c r a , for fermion states 
at site r with spin projection a. The sums in ([!]) run over 
pairs < , > of nearest neighbour sites of the lattice. 

In order to cast the partition function Z = 
irexp{— f3H.} into the form of a path integral, we intro- 
duce the following coherent statesEj for each lattice site 
r (we omit the site index wherever this does not lead to 
ambiguities): 
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Here, 77, 77* are Grassmann variables, and |0> is the vac- 
uum with no spins present. Obviously, |f2 > is a one- 
fcrmion state with spin in the direction (6, </>) on the 
sphere. The inclusion of a phase factor which depends 
on the third Euler angle if) is necessary, if |f2> and 77 are 
to satisfy the following requirements: (i) \Q > is to trans- 
form covariantly under all SU(2) operations, while (ii) 
77 remains invariant under these operations. As will be 
seen shortly, both requirements appear natural in view of 
the physical interpretation of the variables fl — (tp,6,(j)) 
and 77. It is easily checked that the states \uj > include 
the empty and the spin-up and spin-down states of the 
fermion, but exclude the doubly occupied state. In fact, 
the states \lu > are complete in this three-dimensional 
space, 



J did \u>xuj\ 



dn^ I d^drie-i v \uj><lo\ = 



J dQ^J dri*dr/e- 2 "^ {17 77* |0 >< 0| + |fixfi|} 
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where J dfl ■ ■ ■ = J Q v dij) d9 sin 9 j Q w d<f> ■ ■ ■ . The 
Euler angles describe the spin degrees of freedom and 
the Grassmann variables 7/ and 77* a spinless fermionic 
hole. Thus the requirement (ii) is cogent. With these co- 
herent states, the standard steps towards a path integral 
representation of the partition function yield 



Z M = J T> M p] J V M [rf,ri\ 
m . 

H <a-l,r|0-,r> -e~ 2 <^ 



exp 
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with Pm[0] = n r =ir^T,r and T> M [r]*,rj\ 



Y[ T —i r drj* T dr\ T r . Here, sums and products run over 
the M (imaginary) time slices of the Trotter decomposi- 
tion of exp(~(3H) (subscript r) and over the N s lattice 
sites (subscript r). At = (3/M is the width of the time 
slices. The boundary conditions are fio,r = ^M,r and 
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Vm't- The final step in the conversion of Zm 
into a path integral would be to take the continuum limit 



M — > 00 in the time direction. This meets with difficul- 
ties. After the expansion of the second exponential in (^) 
to linear order in At, a factor Y[ T r e^ r>T - rT,T - r remains in 
the integrand which is therefore ill defined in the limit 
M — ► 00. However, as we shall show now, the individual 
summands Z n ~ hole in the decomposition 



crO—hole 
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i-?! — hole 
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<jn—hole 
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of the partition function Zm into 0-, 1-, . . . n-hole contri- 
butions are well defined in the T-continuum limit. Since 
in this work we are interested in the one-hole effective 
Hamiltonian, which can be extracted from z 1-11 " 1 ^, we 
confine our attention to the first two terms of (§). The 
crucial elements in this decomposition are the weight fac- 
tors e ~ 2 Vr'?r.r = \ - 2r]* tr r) TI in the integrand of (Q). 
From (||) it follows that the identity in such a factor 
projects onto configurations with a hole at (t, r) while 
the second term, —2 77* r r\ T r , projects onto configurations 
with an electron at (t, r). Thus, by retaining in the inte- 
grand of (H) only the term J7 T r ( — 277*^77^^) of the expan- 
sion of J| T r (l — 2?7* r 77 T r ) , one obtains Z°~ hole , the par- 
tition function with no hole at all which is, of course, the 
partition function of the Heisenberg Hamiltonian. Sim- 
ilarly, by retaining those terms of the expansion which 
contain for each time slice r just N s — \ factors — 277* r 77 Tr 

one obtains Z M ~ hole . (The number of holes remains con- 
stant in time.) Using the identity 
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which holds between the one-hole projectors on both 
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sides, we can now cast Z\,, hole into the form 



jl — hole 
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The additional factor of | associated with each time slice 
t finds a natural explanation: It compensates for the in- 
tegrations overj-the redundant spin variable at the posi- 
tion of the hole&3. Obviously, the integrand of (|7]) is well 
defined in the T-continuum limit M — > 00. Performing 
this limit, one obtains 



rl — hole 



V[Q] I V[q\rf\ e-fo drC 



(8) 



1 — hole 
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with the classical Lagrangian for one hole 

£ = T,l 1 - r >*AT)vr(T)} <n r ( T )\n r (T)> 

r 

+ $>;('X(t-)- <u{t)\H\u{t) >i- hole . (9) 

r 

Here, \lj(t) >= J] r |w r (r) >, and < cd(r)|W|o;(r) >i-hole 
denotes that part of the expectation value which is bilin- 
ear in the Grassmann fields 77* , 77. 

It should be clear from these considerations that in or- 
der to determine the effective interaction between a pair 
of holes, one will have to analyze the 2- hole contribution 
to Zm- In calculating it from the 1-hole contribution by 
integrating over the spin degrees of freedom, one would 
miss a part of this interaction. 

Evaluating (||) with the states (|^) and the Hamilto- 
nian (Q) of the t — J model on an arbitrary lattice, we 
obtain C = Ckin + £t + Cj, where (we omit total time 
derivatives) 



Ckin — ^ 



VrVr + ^rVr^r 



(1 - VrVr)^ COs(0 r ) 
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and 



Cj = -J 
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The Heisenberg term Cj accounts for the interaction be- 
tween spins S = I (sin 9 cos <f>, sin 9 sin cf>, cos 9) at sites 
where no hole is present. Ct describes the hopping of 
the hole; the hopping energy depends on the state of the 
spins. Finally, Ckin contains the kinetic terms of a spin- 
less hole (rfrj), and a spin |, i.e. <ftcos(9). In addition, 
two terms couple the hole to the angular degrees of free- 
dom. The second cancels the kinetic term of the spin 
at sites where a hole is present. The first has the form 
of a gauge interaction with the field ip. It is tempting 
to use a parameterization in which the factor exp(itp) is 
absorbed in 77, since then, all the ip dependent terms in 
C disappear. However, after such a gauge transforma- 
tion, 7/ is no longer invariant under SU(2) operations, as 
explained above in connection with the choice of the co- 
herent states (|^). Since this has unwanted consequences, 
we do not follow this route. 

In summary, our procedure for obtaining effective 
Hamiltonians consists of the following steps: first we 



represent the partition function of the initial quantum 
Hamiltonian as a discrete-time path integral; then, wc 
identify that part of the partition function that corre- 
sponds to the number of holes we wish to consider; in 
this part, we then perform the continuum limit to obtain 
the classical Lagrangian; finally, we translate back to an 
effective quantum Hamiltonian. 

The Lagrangian £, ( ^0[jl2| ), still represents the full non- 
linear problem of a hole interacting with the spin back- 
ground. In the following, we shall confine ourselves to 
a spin wave expansion around the classical spin ground- 
state of C. In this picture, the hole moves in a spin 
background described by angular fields which deviate lit- 
tle from their groundstate. The spin wave expansions 
have, of course, to be performed separately for the case 
of collinear spin order (square lattice AF) and for the 
case of planar spin order (triangular lattice AF). We shall 
briefly rederive the quantum Hamiltonian, Tin , which is 
well known for the case of collinear spin order. For the 
case of planar spin order, we obtain a significantly dif- 
ferent effective Hamiltonian Ha- To check the validity 
of the approximations that lead to TIa we shall calcu- 
late the one hole spectral properties that follow from Ha 
and compare them with results obtained by numerically 
exact diagonalization of the t — J model for small finite 
lattice cells. 



III. ONE HOLE ON A SQUARE LATTICE 

The two sublattice classical AF order is reproduced by 
assigning the value ir/2 to 9 T and the values 4>^ = ±ir/2 
to 4> r for r from the A or B lattice, respectively. The de- 
viations from the ordered groundstate, x r and p r defined 
by 



V22 



7T 1- 
= - + V2 Pr , 



(13) 



are canonically conjugate harmonic oscillator fields as can 
be seen from Ckin- Now, we expand C up to quadratic 
order in the amplitudes a — (x — ip)/^/2 and a* and 
keep in the quadratic term only the leading zero-hole 
contribution. Then, the spin and the hole degrees of 
freedom decouple in Ckin- We redefine the hole field by 
h = 77 e^^ _ ^ (0) ^ . Then, all the phases in the hopping 
term in the path integral become equal; the field ip disap- 
pears from C and is integrated out. Finally, we translate 
back to operator form and get 



Tin = —t /i r i (al — a T i) + h. 

<r,r'> 
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<r,r'> 



t t 

a r , — a [ r a v , — a r ,a r 
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Tin is, up to a canonical transformation (a— — > — ia r ) 
identical with the Hamiltonian considered inci. The last 
term in 7ip , containing the hole number, accounts for the 
breaking of four bonds per hole in the classical ground- 
state. 



and Sj — (cos ip j, sin ip j), ipj — ^-j, are the vectors to 
the three nearest neighbour sites on the triangular lat- 
tice. As opposed to the square lattice, the hole can hop 
on the triangular lattice without emission or absorption 
of a magnon. This is expressed by the second line of jl5| ) 
where 



IV. ONE HOLE IN THE TRIANGULAR LATTICE 

In the classical limit, the groundstate of the t—J model 
at half filling is the well known planar 120° spin struc- 
ture: ei Q) = tt/2 for all r and c^ 0) = for the 
A, B and C sublattice, respectively. As in the case of 
the square lattice, we describe the deviations from this 
ordered groundstate by the spin wave amplitudes x T and 
p r which are canonically conjugate harmonic oscillator 
fields. 

Proceeding as in the square lattice case, i.e. expanding 
to second order in the amplitudes a and a* and integrat- 
ing over ijj, one would be left with terms in £ t that couple 
the hole motion to zero energy spin modes. But these 
modes are uniform spin rotations which are exact sym- 
metries of the initial model ([!]). Therefore, they cannot 
appear in any order of a perturbation expansion. In the 
Appendix, we trace this apparent inconsistency back to 
the integration over the field ip. We show how to derive 
a result which is free of zero modes, as it should be, and 
we explain why the same difficulty does not occur in the 
case of the square lattice. Retaining terms of the same 
order in the amplitudes a, a* as in the square lattice case, 
we find the following effective Hamiltonian for one hole 
in the triangular AF (in momentum space) 
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The terms in the first line are nothing but the linear spin 
wave approximation (LSW) of the triangular Heisenberg 
AF. The last line contains the coupling between the hole 
and the magnons with the hole-magnon vertex 



jq,q' = 3V3i 



J 

-*7q' + 7^q-q' 



^3 (^+Wq-q' 



*7q-c 



t(Wq - Wq>) + -(1 - Wq-q') 



(16) 



Here, Wq is defined by 



3J 

T 



3tw a 



(18) 



is just half the dispersion of a particle hopping between 
nearest neighbour sites of the empty triangular lattice. 
Ti/x constitutes our main result for the Hamiltonian of 
one hole in the t — J model on a triangular lattice in 
spin wave approximation. In the following, we analyze 
its content by treating the hole magnon coupling in the 
simplest approximation. 



A. Born approximation of the selfenergy 




FIG. 1. Spectral density A(k,u) in SCBA approximation 
for N s = 36 x 36 sites, t = J = 1, and for k values taken along 
a straight line from F to K (cf. Fig. §) (tj = 0.1). 

Following previous treatments of the square lattice 
caseo, we employ the SCBA to obtain the selfenergy 
E(k, ui) of one hole. We denote the Green's function of 
one hole by G(k,w) = (w - e k - S(k, w))" 1 . Then, 
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SCBAnJ 



(k ,u- Wk-k' 



(19) 



»q + «7q = 



(17) Here, w q = + 2wq) (1 — Wq) is the spin wave en- 

ergy. Uq and Vq are the Bogoliubov amplitudes, 
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From the Green's function, the spectral density is ob- 
tained as A(k,uj) = — - lim ?? ^o+ LmG(k, u + irj). The 
parameter rj regularizes the spectral density. 
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-5.25000 


0.50000 


0.40000 


4.00000 


-5.62500 


0.66884 


0.31627 


3.18974 
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-7.32396 


1.15197 


0.38544 


4.09035 


-7.25658 


0.64520 


0.32854 


3.39377 
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-11.78091 


0.56862 


0.39986 


4.06326 


-11.82919 


0.65145 


0.32504 


3.35476 


24 1 " 


-12.93870 


0.92897 


0.29581 


3.42937 


-13.04174 


0.65227 


0.32561 


3.36573 


27 


-15.12597 


0.59201 


0.38805 


3.99864 


-15.04634 


0.65052 


0.32572 


3.36510 



TABLE I. Exact and approximate groundstate energies of samples of N s sites of the t — J model with rih(~ 0, 1) holes. Eq: 
groundstate energies in the half filled case. Ei t i oc : energies for creating a localized hole. 5E\: energy gain of the hole due to 
derealization. 



( |i9| ) has been solved by numerical iteration starting 
with E = for a lattice of N s — 36 x 36 sites with peri- 
odic boundary conditions. In Fig. [l] we display for t = 1, 
J = 1 the spectral density A(k, uj) on k-points along a 
straight line from the centre T to the corner K of the 
Brillouin zone (cf. Fig. ||). The main feature of the spec- 
tra is a pronounced quasiparticle peak. For t — 0, the 
quasiparticle dispersion is flat throughout the Brillouin 
zone. For increasing t, a dispersion with the minimum at 
r emerges and the peak broadens. Around T, the quasi- 
particle peak persists for values of t up to 10 J; near K it 
already decays for t > 0.2 J. 

B. Comparison with exact numerical diagonalization 
results 

As we have pointed out in the course of the derivation, 
the effective Hamiltonian 7Ya is an approximate result, 
and the computation of the one-hole Green's function 
is based on the self-consistent Born approximation. In 
order to have a check on the validity of these approxi- 
mations, we computed the spectral density A(k, uj) and 
related quantities directly for the t — J model on small 
periodic samples of the triangular lattice with N s sites. 
Numerically exact results for the quantities in question 
wer£ obtained by applying the Lanczos technique (see 
e.g.EJ) to the Hamiltonian of these samples. In these com- 
putations, the exchange constant J has been set equal 



to unity throughout. Results for the groundstate ener- 
gies for half filling (rih — 0) and for one hole (n^ = 1) 
are displayed in Table | together with the approximate 
LSW and SCBA energies for the same sample sizes. For 
t =/= 0, the groundstate wave function of the hole is an 
extended state with crystal momentum at the high sym- 
metry points of the Brillouin zone; T for odd N s and K 
for even N s . 




FIG. 2. The Brillouin zone of the N s — 21 system shows 



a sixfold symmetry. 



^For N 3 = 24 the crystal momentum of the one hole ground- 
state is an interior point close to K of the Brillouin zone. This 
is the reason, why for N s = 24 the values of 8Ef x deviate sig- 
nificantly from the rest. 
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For the half filled case, the t — J Hamiltonian reduces 
to the Heisenberg Hamiltonian, and it is well known that 
in this case the LSW results Eq SW are in good agree- 
ment with the exact results Eq x , in particular for the 
larger system sizestll. Turning to the one-hole results, we 
first discuss the case t — in which the hole is localized. 
Exjac is the energy needed to create a localized hole in 
the half filled groundstate. There is a strong even-odd 
staggering in the exact results (2nd column of Table |) : 
in the even samples (N s even) in which the spins pair 



up to a singlet, S to t = 0, the 'binding energy' per elec- 
tron is larger than in the odd systems which contain an 
unsaturated spin, Stot — 1/2. This is a finite size effect 
which, as is seen in Table Q, decreases very slowly with the 
sample size. Nevertheless, the exact values Ef x loc appear 
to converge to the SCBA value EffJ? A in an alternat- 
ing fashion. That the even-odd staggering is absent from 
Eif c A i s understandable, since the SCBA is based on 
the spin wave spectrum of the infinite sample. 
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FIG. 3. Spectral density a) at the V point and b) at the corner K of the Brillouin zone for t = 0.1 and t = 1 for J = 1. 
The solid line shows the exact results for the N s = 21 sites, the dashed (dotted) line the result of the SCBA for N s = 21 
(N s = 36 x 36) sites, uo is measured against Eq x and Eg sw (rj = 0.1). 



For finite hopping amplitude t ^ 0, the hole becomes 
delocalized. ft is suggestive to write the hole energy as 
Ei = Eijoc — 8E\(t), where 6Ei(t) is the energy gain 
due to delocalization of the hole wave function. The re- 



sults in Table [j] show that 8E\ increases linearly with t, 
SEi(t) = at. This assertion has been verified by com- 
puting SEi(t) for a sequence of values 0.1 < t < 1 not 
shown in Table B The coefficient a is seen to differ signif- 



G 



icantly between the exact results and the SCBA, a ex « 4, 
a sc'BA ^ 3 35 At present, we are not in a position to 
decide whether this discrepancy is inherent in the effec- 
tive Hamiltonian or whether it is a deficiency of the 
SCBA. 

Next we compare the exact and the SCBA results for 
the spectral density j4(k, u>). Fig. || shows the spectral 
density of one hole in the N s = 21 sample at the cen- 
tre (r) and at the corner (K) of the Brillouin zone (cf. 
Fig. H) for the parameter values J — 1 and t = 0.1, 
t = 1. The solid lines are the results of the exact diago- 
nalizatiori|-which were obtained from a continued fraction 
expansioiJla of the Lanczos results. The dashed lines are 
the result of the SCBA for N s = 21. For t = 0.1 there is a 
well defined quasiparticle peak in the exact results at the 
bottom of the spectrum. This peak is nicely reproduced 
by the SCBA. 

For t= 1, a quasiparticle peak is still visible at T while 
one finds a broad structure at the corner of the Brillouin 
zone. The SCBA reproduces the peak in the spectrum at 
T; however, its position is shifted to higher energies. At 
K, the SCBA spectrum shows a structure which becomes 
broader when the system size is increased from N s = 21 
to N s — 36 x 36 (dotted lines). On the other hand, in 
all cases where we find a quasiparticle peak, the SCBA 
results for N s — 21 and for N s — 36 x 36 are indistinguish- 
able. Thus we conclude that for t = 1 the quasiparticle 
peak disappears if one changes the wavevector from the 
r point to the corner of the Brillouin zone, cf. Fig. [I]. 
While there remain differences, we find that overall the 
SCBA reproduces the qualitative features of the exact 
results. 



We define the dispersion as the difference of the position 
of the maximum of A(k, uS) for a given k and the one 
hole creation energy E^f a ^ A (see Table |), i.e. the po- 
sition of the maximum for t = 0. In Fig. ^, we show 
the dispersion for AT S = 21 sites for t = 0.1 and J = 1. 
The circles (connected by solid lines) show the results of 
the exact diagonalization. They follow closely the free 
dispersion of the hole, e q — 3 J/4 = — 3tu> q (shown as 
boxes). The SCBA (diamonds) yields an improvement 
towards the exact results for all points of the Brillouin 
zone. Obviously, the free dispersion is already a good 
approximation of the exact results. 



V. SUMMARY 



In this work, we have developed a formalism which al- 
lows the derivation of effective Hamiltonians for the mo- 
tion of holes in an arbitrarily ordered spin background. 
Fundamental to this development was the construction 
of a coherent state approximation for the electron states, 
which automatically excludes double occupancy. In this 
representation, the partition function of an arbitrary 
number of holes can be cast into the form of a path inte- 
gral which still contains the full non-linear coupling of the 
holes to the spin degrees of freedom of the background. 
As a first application of this formalism, we have rederived 
the well known effective Hamiltonian H a for the motion 
of a single hole in the collinear antiferromagnetic order 
of the square lattice antiferromagnet in the linear spin 
wave approximation. The derivation of the effective 
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FIG. 4. The dispersion of one hole as denned in the text 
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Finally we discuss the energy dispersion of one hole. 
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Hamiltonian Ti/s of a single hole moving in the planar 
spin arrangement of the triangular antiferromagnet also 
makes use of the linear spin wave approximation. How- 
ever, in this case, special care has to be taken to ensure 
a proper treatment of the Goldstone modes which must 
not couple to the hole. In contrast to the square lattice 
AF, where hopping of the hole from a lattice site r to a 
neighbouring site r-\-8 is necessarily accompanied by a 
spin flip at site r, the same hopping process requires only 
a rotation of the spin at r by ±120° in the case of the 
triangular AF, see Fig |^. 

In other words, if l^r "* > is the groundstate with one 
hole at r in the undistorted spin background of the AFs, 
then the hopping matrix element <^f^ s \c r +saC 
is always zero for the square lattice AF, while it is +5 for 
the triangular AF. This accounts for the main difference 
between the hole dynamics in the collinear square lattice 
AF and the planar triangular AF. 

On the square lattice, the hole can hop due to the 
zero point fluctuations of the AF background. In this 
case, magnon-assisted hopping happens predominantly 



between the sites of the same magnetic sublattice. This 
leads to a minimum of the hole dispersion at a quarter 
of a reciprocal lattice vector of the square lattice. By 
contrast, on the triangular lattice, hole hopping between 
nearest neighbour sites, i.e. between different magnetic 
sublattices, is also possible without magnon assistance. 
Thus, for J >> t, where magnon assisted processes are 
energetically suppressed, the hole dispersion in the tri- 
angular AF will be dominated by the bare (unassisted) 
hopping processes. This has been confirmed by compar- 
ing the exact dispersion of the bare hole in the 21 site 
sample for t = 0.1J. The approximate inclusion of the 
magnon-assisted processes in the calculation of the dis- 
persion by the self-consistent Born approximation has 
been found to improve the agreement with the exact dis- 
persion. The exact one-hole Lagrangian jlO|-p^[), eval- 
uated in spin wave approximation, leads to results in 
reasonable agreement with exact numerical diagonaliza- 
tion of small systems. C provides thus a firm basis for 
a derivation of a continuum model which is capable of 
describing the full non linearities of the spin fields in the 
presence of single holes. 



APPENDIX A: 

A straightforward expansion of C around the classically 
ordered state with respect to the fields 9 and (p, and a 
subsequent integration of ip leads to the appearance of 
zero energy spin wave amplitudes in linear order in the 
hopping term C t . This is in contradiction to the spin ro- 
tational invariance of the model (§). In order to resolve 
the contradiction, we now verify that this invariance per- 
sists in an appropriate spin wave expansion. In the first 
order expansion around the state given by (p r = (pi , 
r = § (cf • |l^) , the critical term Ct reads 



|Mv-V>r) 
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.Pr Pr' 
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This expression is valid for both the square, and the tri- 
angular lattice. We now perform an infinitesimal homo- 
geneous rotation in spin space parametrized by the angles 
e. The corresponding transformation of the fields can be 
read off from (||) . The field 77 is unchanged and the linear 
changes in ip, 6, and <f> are (we express the last two fields 
by x and p) 



< = ^ r + 2 
x' r = x T + v2e 
p' = Pr + V2 



e cos 



do) 



-e sin 



e v sm«!>i 0) 



+ € y COS 



(A2) 



Substituting these rotated fields for the original ones in 
(Al), one sees that e z drops out and the change in the 
field p r induced by the other two rotations e x and e y is 
compensated by the corresponding change in the field ip r - 
Thus, we see that Ct, and hence the whole Lagrangian, 
is invariant under an infinitesimal rotation in spin space 
as it should be. It is only after the integration over the 
field ip that the invariance appears to be lost. 

Turning to the case of the square lattice, where (pi^ = 
±tt/2 (cos4 0) = 0), the zero modes are seen to decouple 
in ip T and p r (cf. (A2)), and in Ct the term of zeroth 
order in the fields a; r and p r vanishes. It is easily verified 
from the spin wave expansion of Cj that the zero modes 
are given by x r — x, p r = and x r = 0, p r — psin^ ^ 
which both cancel in Ct- After substituting 77 e^^ - * ° ) 
by h, we can integrate over tp, and return to the operator 
form. Thus, the spin wave expansiop_.as sketched in the 
main text leads to the correct resulto. 

In the case of the triangular lattice, where ^r "* = 
3p 0, — ^ for the three sublattices respectively, the zero 
modes in ip r and p r mix. Again, it is easily verified from 
the spin wave expansion of Cj that they are given by 
x r = x, p r = 0, and x r = 0, p r = psin</>r°\ and x r = 0, 
p r = p cos 4^ ■ The first mode cancels in Ct , but the last 
two remain, if one disregards ip. However, we are free to 
apply a canonical transformation to the fields ip, x, and 
p in the path integral. We choose 
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X r -> X T + \ - 7] r T], 



Pr -> Pr 



3=1 



(A3) 



(<Jj = (cos^j, siniy9j), tpj = -j-j are the three nearest 
neighbour vectors.) With this choice (i) the new field ip 
remains invariant under homogeneous rotations in spin 



space, i.e. it does not contain a zero mode, and (ii) Ckin 
is changed only by a total time derivative which yields 
no change in the action (the latter condition is equivalent 
t o t he transformation being canonical). After applying 
([A3]) to C, all zero modes cancel in £ t . We substitute 



■q e 



■ (i/> 4> ( ) b v h, integrate over tp, and return to the op- 



erator representation. The result for Tl/\ is the expression 
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